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We report the storage of microwave pulses at the single-photon level in a spin-ensemble memory 
consisting of 10^° NV centers in a diamond crystal coupled to a superconducting LC resonator. 

The energy of the signal, retrieved 100 jis later by spin-echo techniques, reaches 0.3% of the energy 
absorbed by the spins, and this storage efficiency is quantitatively accounted for by simulations. 

This hgure of merit is sufficient to envision first implementations of a quantum memory for super¬ 
conducting qubits. 


Superconducting qubits are attractive candidates for 
solid-state implementations of quantum information pro¬ 
cessing, but suffer from coherence times shorter than 
^ 100/rs [1-3]. To circumvent this issue, it has been pro¬ 
posed to use ensembles of spins in semiconductors [4-9] 
as a multimode quantum memory, able to store multi¬ 
ple qubit states over longer periods of time, and to re¬ 
trieve them on-demand [10]. Inspired by research on op¬ 
tical quantum memories [11-13], realistic protocols have 
been proposed recently [14, 15]. The state of a super¬ 
conducting qubit is first converted into the state of a mi¬ 
crowave photon, propagating or trapped in a resonator. 
This photon is then resonantly and collectively absorbed 
by the spin ensemble, resulting in a transverse magneti¬ 
zation which, due to the spread of resonance frequency 
within the ensemble, decays in a time T 2 called the free- 
induction decay (FID) time. The write step is later fol¬ 
lowed by the need to read the stored quantum state. Both 
protocols [14, 15] propose to apply sequences of tt pulses 
to the spins, combined with dynamical tuning of the res¬ 
onator frequency [16, 17] and quality factor [18, 19] in 
order to trigger the rephasing of the spins, resulting in 
the emission of an echo at a chosen time that faithfully 
reproduces the initial quantum state. 

Whereas the transfer of a qubit state into a spin en¬ 
semble has been demonstrated experimentally [20-23], 
implementing the read step remains the major obstacle 
to an operational microwave quantum memory. An in¬ 
termediate goal consists in storing a classical microwave 
pulse with an ultra-low power corresponding to an aver¬ 
age of 1 photon in the resonator and to retrieve it as an 
echo after a refocusing pulse, as was achieved at optical 
frequencies [24, 25]. First results in this direction were 


obtained using ensembles of negatively-charged nitrogen- 
vacancy (NV) colour centers in diamond [26, 27] and of 
rare-earth ions in a ¥28105 crystal [28]. The NV’s elec¬ 
tronic spin is a spin triplet {S = 1 ) well suited for a 
quantum memory because of its long coherence times in 
pure crystals [29, 30] and the possibility of repumping it 
into its ground state = 0 by optical irradiation at 
532nm [31] (see Fig. 1). In [27], successive low-power 
microwave pulses were stored in an NV ensemble, and 
retrieved later as a series of echoes after a refocusing 
microwave pulse was applied. A key aspect in this ex¬ 
periment was an active reset of the NVs to increase the 
repetition rate of successive experimental sequences to 
obtain sufficient statistics; this was achieved by applying 
optical pumping laser pulses injected through an optical 
fiber introduced in the cryostat. The echo efficiency, de¬ 
fined as the ratio of the echo energy and the stored pulse 
energy, was, however, not sufficient in [27] to observe an 
echo below 100 photons on average in the resonator. 

Here, using a sample with a longer coherence time and 
an improved optical pumping scheme, we increase the 
echo efficiency and the storage time by one order of mag¬ 
nitude. This allows us to observe an echo with an initial 
pulse power corresponding to on average only one photon 
in the resonator. The diamond single crystal was synthe¬ 
sized by the temperature gradient method at high pres¬ 
sure and high temperature (HPHT) using 99.97% ^^C- 
enriched pyrolytic carbon prepared from ^^C-enriched 
methane as a carbon source [32], resulting in a nomi¬ 
nal 300 ppm concentration of nuclei. The original 
1.4 ppm concentration of substitutional nitrogen impuri¬ 
ties (so-called PI centers) was partially converted into 
NV centers by 2MeV electron irradiation at room tern- 
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FIG. 1. (a) Experimental setup and principle of the exper¬ 
iment. An ensemble of ~ 10^° spins is inductively coupled 
to a planar superconducting LC resonator of frequency ujr 
(with a collective coupling constant gens)-, cooled at 10mK. 
The resonator is measured in reflection through an input cou¬ 
pling capacitance. Microwave pulses are produced by mix¬ 
ing a continuous microwave source with dc pulses generated 
by an arbitrary waveform generator (AWG). They drive the 
spins via the microwave current induced in the resonator in¬ 
ductance. The reflected microwave signal (including the emit¬ 
ted echo) is amplified at low-temperature and demodulated 
at room-temperature, yielding its amplitude A(t) and phase 
0(t). (b) The spins are Nitrogen-Vacancy color centers in dia¬ 
mond, which consist of a nitrogen impurity next to a vacancy 
of the diamond lattice. A dc magnetic field applied par¬ 
allel to the chip along the [110] crystalline axis so that only 
NV centers whose axis are non-orthogonal to the field (shown 
in blue on the figure) are Zeeman-shifted and contribute to 
the signal. Laser pulses can be sent onto the diamond via a 
direct optical access to the cryostat mixing chamber, (c) NV 
centers energy levels in a weak magnetic field. The electronic 
ground state is a spin triplet iS = 1, with a zero-field split¬ 
ting Djl'K — 2.88 GHz, coupled by hyperfine interaction to 
the nuclear spin triplet 7=1. This splits each of the 
= 0) ^ \ms — +1) transitions into a triplet of lines. 



FIG. 2. (a) Phase 0 of the resonator reflection coefficient 

measured with a vector network analyzer (blue dots), 
yielding — 2.915 GHz and Q — 650 (red line is a fit to 

the data), (b) *S'ii(ct;) as a function of Bjyy around 1.8 mT, 
showing the NV centers as a triplet of absorption dips. The 
inset shows IaShKHatv) at Ct;/ 27 r = 2.915 GHz; blue dots are 
data, and red line is a fit to a sum of three Lorentzians with 
linewidth 0.012 mT. (c) Measured spin polarization for a laser 
pulse of power Pl = 280/iW as a function of its duration Tl, 
renormalized to its maximal value. Dashed red and black lines 
indicate the laser pulse durations Tl — 0.2 and 0.1s used in 
the experiments shown in Figs. 3 and 4, corresponding to 
relative polarizations pipmax = 0.72 and 0.62. 


perature followed by annealing for 2 hours at 1000° C, 
yielding a final [NV~] = 0.4ppm and [PI] = 0.6ppm. 
A scheme of the experimental setup is shown in Fig. la. 
The diamond is glued onto the inductance of a supercon¬ 
ducting planar lumped-element LC resonator patterned 
in a niobium thin-film on a silicon substrate. Microwave 
pulses can be sent to the resonator input; the reflected 
signal is amplified at 4 K and its amplitude and phase are 
measured at room-temperature by a homodyne detection 
setup. For resetting the spins, laser pulses at 532 nm can 
be sent onto the sample through direct optical access in 
the dilution cryostat. 

The resonator reflection coefficient 5 'ii(cj), measured 
with a network analyzer, is shown in Fig. 2a, yielding 
the resonance frequency Url^ir = 2.915 GHz and quality 


factor Q = 650, fixed by the coupling to the measure¬ 
ment line through the input capacitor. NV centers are 
detected by their absorption of the microwave whenever 
their transition frequency matches uJr- The energy levels 
of the NV centers are schematically shown in Fig. Ic. The 
electronic spin is coupled by hyperfine interaction to the 
spin-triplet (7 = 1) nuclear spin of the atom, result¬ 
ing in a splitting of the \ms = 0) ^ \ms = +1) transition 
into three resonances separated by 2.2 MHz correspond¬ 
ing to the three different mj states of the A dc 

magnetic field B^v is applied parallel to the chip, along 
the [110] crystalline axis of the diamond. Out of the four 
possible orientations of NV centers along [111] crystalline 
axes, two are perpendicular to the field and are there- 
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fore not Zeeman-shifted, so that they do not contribute 
to the signal. The remaining two families are brought 
into resonance with at ^ 1.8 mT, as shown in 

Fig. 2b where the hyperfine triplet is clearly seen as dips 
in |5'ii(cj)| when they cross the resonance. A ^ 200 kHz 
Full-Width Half Maximum (FWHM) linewidth is mea¬ 
sured for each line of the triplet, much narrower than in 
previous work [20, 27], due to a lower PI center concen¬ 
tration and to the isotopic enrichment in ^‘^C. 

For optical pumping of the NV centers, the laser beam 
is focused to a 0.6 mm diameter at the sample level, with 
a power of 0.28 mW. In this geometry, it is straightfor¬ 
ward to optimize the laser beam position on the sam¬ 
ple in order to minimize the amount of power needed 
to reset the spins in their ground state. The efficiency 
of the optical pumping is measured as explained in [27]. 
The experimental sequence includes an initial strong mi¬ 
crowave pulse that saturates the spins, followed by a 
laser pulse of varying duration T^. After a 300/is delay 
necessary for relaxation of the quasiparticles generated 
in the superconducting thin film and in the silicon sub¬ 
strate, the reflected amplitude of a few-photon microwave 
pulse reveals the spin polarization. The extracted po¬ 
larization level p{Tl) is shown in Fig. 2c. It changes 
only slightly above Tl = Is, indicating that the maxi¬ 
mum NV polarization possible with irradiation at 532 nm 
{Pmax = 90% according to [33]) is reached. Compared 
to earlier work [27] where the laser position could not 
be optimized, the maximum polarization can be reached 
with ^ 20 times lower pulse energy. This makes it possi¬ 
ble to perform the experiments at a faster repetition rate 
(0.2 Hz) and at a cryostat temperature of 100 mK instead 
of 400mK [27]. 

High-power Hahn echoes are measured at B^y = 
1.74 mT, using microwave pulses ai uj = ujr according 
to the sequence shown in Fig. 3a. The sequence starts 
with a laser pulse of duration Tl = 0.2 s, resulting in a 
spin polarization p = 0.72 pmax = 0.65. At t = 0 , a first 
pulse generates a transverse magnetization in the ensem¬ 
ble, followed by a refocusing microwave pulse at t = r 
which induces rephasing of the spins at 2 r and emission 
of a spin-echo into the measuring line, as seen in Fig. 3a. 
Note that due to spatial inhomogeneity of the microwave 
field generated by the planar inductance, it is not possible 
to apply a well-defined Rabi angle to all the spins, which 
results in a reduced echo visibility. The echo amplitude 
is measured as a function of the delay 2 r between the 
first pulse and the echo, and is found to decay approx¬ 
imately exponentially with a time constant T 2 = 84/is 
(see Fig. 3b). Decoherence occurs due to dipolar inter¬ 
actions with the bath of paramagnetic species present 
in the sample nuclei, PI centers, and NV centers), 
whose dynamical evolution causes a randomization of the 
phase acquired by NV centers during the two halves of 
the spin-echo sequence. The nuclei bath precesses at 
the Larmor frequency '^uBnv = 27r'I30kHz ( 7 ^ being 
the gyromagnetic ratio), giving rise to a characteris¬ 
tic oscillation pattern [34-36] in the spin-echo amplitude. 
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FIG. 3. (a) Hahn echo sequence. The spins are first reset 
in their ground state by a laser pulse of power 280/iW and 
duration 0.2 s. A first Ifis microwave pulse ^ at Ct;r, of power 
—71 dBm, induces a transverse magnetization which decays 
within T 2 . A I/is-long microwave refocusing pulse (R) of 
power —20 dBm is applied at r = 50//S, which rephases the 
spins at 2t. The microwave amplitude (blue curve) shows 
both the reflection of the two microwave pulses driving the 
spins (with their amplitude trimmed by saturation of our de¬ 
tection chain, as indicated in red), as well as the echo emit¬ 
ted at 2t upon rephasing of the spins, (b) Measured decay 
of the echo amplitude A as a function of 2r (open circles). 
Calculated decay due to a bath of 213 ppm of (dashed or¬ 
ange curve), 0.6ppm of PI centers causing spectral diffusion 
(dashed green line), 0.2ppm of NV centers causing instanta¬ 
neous diffusion (dashed blue line), and the combination of the 
three contributions (solid red line). The theory curves have 
been scaled in amplitude according to the data. An exponen¬ 
tial fit (black solid line) yields a coherence time T 2 = 84/is. 


visible in the data of Fig. 3b. The dynamics due to flip- 
flop events within the PI center bath is responsible for a 
decoherence process knwon as spectral diffusion [37]. Fi¬ 
nally, the bath consisting of NV centers at frequency ujr 
(only half of the total NV concentration) unavoidably un¬ 
dergoes spin flips due to the refocusing pulse itself, which 
constitutes an efficient decoherence process called instan¬ 
taneous diffusion [38]. The various contributions of each 
bath were calculated using the cluster-correlation expan¬ 
sion method [39, 40], with concentrations [PI] = 0.6 ppm, 
[NV~] = 0.4/2 = 0.2 ppm, and [^^C] = 213 ppm, com¬ 
patible with the sample parameters. Good agreement 
with the data is obtained (see Fig. 3b). Overall, dipolar 
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FIG. 4. Spin-echo at the few-photon level. The experimental 
sequence is the same as in Fig. 3, but with an initial microwave 
pulse having a power corresponding to (a) 60 and (b) 1 photon 
inside the resonator. The signal was averaged over 3 x 10^ (a) 
and 5 x 10^ (b) sequences, with a repetition rate of 5 Hz (a) 
and 10 Hz (b), limited by the laser pulse duration. Blue solid 
lines are experimental data, red dashed-dotted lines are the 
results of simulations as explained in the text. 


interactions between NV centers appear to be the domi¬ 
nant source of decoherence in our experiment. 

Since the echo efficiency was limited by the finite spin 
coherence time in earlier work [27], a significant improve¬ 
ment is expected with this new sample. The echo effi¬ 
ciency is first measured by performing a Hahn echo se¬ 
quence with a low-power microwave pulse. The experi¬ 
mental sequence, shown in Fig. 4, starts with a 0.2 s re¬ 
pumping laser pulse, followed after 300 /rs by a microwave 
pulse populating the resonator with on average 60 pho¬ 
tons, and, r = 50/is later, by a refocusing pulse. A spin- 
echo is detected at t = 2r. The efficiency, defined as the 
energy recovered during the echo divided by the absorbed 
energy, reaches = 0.3%. 

The whole spin-echo sequence was simulated by numer¬ 
ical integration of the coupled differential equations ob¬ 
tained after discretizing the resonance frequency and cou¬ 


pling constant distribution in the spin ensemble (see [14]). 
A 130 kHz linewidth of each HF peak and a collective 
coupling constant ^ens/^Tr = 410 kHz (when all spins are 
polarized), are compatible with the experimentally deter¬ 
mined parameters and yield quantitative agreement with 
the shape of both the absorbed microwave pulse and the 
spin-echo amplitude. The discrepancy noted in [27] is 
absent in the present experiment, probably because de¬ 
coherence is negligible during the driven evolution. The 
finite T 2 and the imperfect tt pulse due to the spread in 
Rabi frequencies are the main factors limiting while 
the finite cooperativity C = 0.22 [14] limits both the ab¬ 
sorption and the echo emission out of the cavity. 

Owing to the larger value of F^, it becomes possible to 
reach the level where a spin-echo can be observed for an 
initial pulse populating the resonator with only a single 
microwave photon on average. This is shown in Fig. 4b. 
Note that a shorter repumping time of 0.1s (with the 
same laser power) was used, in order to enable a larger 
number 5 x 10^ of repetitions of the experiment. The 
shorter repumping step yields a lower spin polarization 
p = 0.56 as shown in Fig. 2 and a lower cooperativity 
of 0.19, which results in a correspondingly lower echo 
efficiency than in Fig. 4a. These results are again quan¬ 
titatively reproduced by the simulations with the same 
parameters as mentioned earlier, using the experimen¬ 
tally determined repumping efficiency. 

The coherence times demonstrated in this experiment 
match those requested in a realistic quantum memory 
protocol [14], which suggests that a first implementation 
is within reach. The remaining challenges are the im¬ 
provement of the refocusing pulse using adiabatic pas¬ 
sage as demonstrated recently [41], and the integration 
of dynamical tuning of the resonator frequency and qual¬ 
ity factor with more elaborate spin-echo sequences. The 
latter is needed in particular to silence the echo emis¬ 
sion [11, 12] in-between the two tt pulses in the course 
of the read step of the procotol. The microwave currents 
needed to drive the spins during the refocusing pulses are 
much stronger than typical Josephson junctions critical 
currents. This precludes the use of integrated SQUIDs as 
in [20, 42], while a combination of coupled linear and tun¬ 
able resonators [43] may be employed as tuning elements 
in the resonator. 

In conclusion we report the measurement of a spin-echo 
with an initial microwave pulse at the single-photon level. 
The figures of merit reached are sufficient to envision first 
implementations of a spin-ensemble multi-mode quantum 
memory for superconducting qubits. 
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